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LONG-TERM  GOALS 

Couple  together  analytical  and  numerical  modeling  of  oceanographic  and  surface  wave  processes, 
acoustic  propagation  modeling,  statistical  descriptions  of  the  waveguide  impulse  response  between 
multiple  sources  and  receivers,  and  the  design  and  perfonnance  characterization  of  underwater 
acoustic  digital  data  communication  systems  in  shallow  water. 

OBJECTIVES 

Develop  analytical/numerical  models,  validated  with  experimental  data,  that  relate  short-tenn 
oceanographic  variability  and  source/receiver  motion  to  fluctuations  in  the  waveguide  acoustic 
impulse  response  between  multiple  sources  and  receivers  along  with  new  communication  receiver 
algorithms  that  exploit  these  channel  characterizations. 

APPROACH 

The  focus  of  this  research  is  on  how  to  incorporate  an  understanding  of  short-term  variability  in  the 
oceanographic  environment  and  source/receiver  motion  into  the  design  and  performance 
characterization  of  underwater  acoustic,  diversity-exploiting,  digital  data  communication  systems.  The 
underlying  physics  must  relate  the  impact  of  a  fluctuating  oceanographic  environment  and 
source/receiver  motion  to  fluctuations  in  the  waveguide  acoustic  impulse  response  between  multiple 
sources  and  receivers  and  ultimately  to  the  design  and  performance  characterization  of  underwater 
acoustic  digital  data  communication  systems  in  shallow  water. 

The  two  major  thrusts  of  this  work  include  participation  in  a  shallow  water  acoustic  communication 
experiment  in  June -July  2011  along  with  subsequent  analysis  of  the  experiment  data. 

KAMI  1  Experiment  (2011) 

A  shallow  water  acoustic  communications  experiment  (KAMI  1)  was  carried  out  in  June- July  2011  off 
the  western  side  of  Kauai,  Hawaii  at  the  Pacific  Missile  Range  Facility  (PMRF).  Both  fixed  and  towed 
source  transmissions  were  carried  out  to  multiple  receiving  arrays  over  ranges  of  approximately  1-14 
km.  Substantial  environmental  data  was  collected  including  water  column  sound  speed  structure 


1 


(CTDs  and  thermistor  strings),  sea  surface  directional  wave  field  (waverider  buoy),  and  local  wind 
speed  and  direction.  The  focus  was  on  fluctuations  over  scales  of  a  few  seconds  to  a  few  tens  of 
seconds  that  directly  affect  reception  of  a  data  packet  and  packet-to-packet  variability. 

Algorithm  Design  and  Experiment  Data  Analysis 

Communication  receiver  algorithm  design  for  shallow  water  is  challenging  due  to  the  likelihood  of 
encountering  highly  dynamic  environments.  In  fact,  one  of  the  reasons  for  the  KAMI  1  site  selection 
was  the  variability  of  the  environmental  conditions  (both  sea  surface  and  water  column)  from  fairly 
calm  to  highly  fluctuating.  One  emphasis  of  our  work  is  mitigating  Doppler  and  Doppler  spread  (both 
due  to  source-receiver  motion  as  well  as  a  fluctuating  sea  surface).  Another  emphasis  of  our  work  is 
multi-user  receiver  design  in  the  presence  of  unsynchronized  users  (to  accommodate  networked 
communication).  Lastly,  a  third  emphasis  of  our  work  is  communication  algorithm  design  for  very  low 
signal-to-noise  ratio  (SNR)  operations. 

WORK  COMPLETED 

A  shallow  water  acoustic  communications  experiment  (KAMI  1)  was  conducted  in  early  summer  2011 
off  the  western  side  of  Kauai,  Hawaii  [8].  Both  fixed  and  towed  source  transmissions  were  carried  out 
to  multiple  receiving  arrays  over  ranges  of  1-8  km  along  with  additional  towed  source  transmissions 
out  to  14  km  range.  The  acoustic  transmissions  were  in  three  bands  covering  3.5  to  35  kHz.  Substantial 
environmental  data  was  collected  including  water  column  sound  speed  structure  (CTDs  and  thermistor 
strings),  sea  surface  directional  wave  field  (waverider  buoy),  and  local  wind  speed  and  direction.  The 
focus  was  on  fluctuations  over  scales  of  a  few  seconds  to  a  few  tens  of  seconds  that  directly  affect  the 
reception  of  a  data  packet  and  packet-to-packet  variability.  The  experiment  region  exhibited 
substantial  daily  oceanographic  variability. 

Analysis  of  the  KAMI  1  experiment  data  this  past  year  has  focused  on  fixed  source  transmissions. 
Specifically,  data  investigating  the  influence  of  surface  interacting  ray  paths  was  analyzed. 

Publications  related  to  this  research  involving  analysis  of  KAMI  1  (and  previous  KAM08)  data  include 
[1-18]. 

RESULTS 

During  KAMI  1,  two  fixed  source  and  vertical  receive  array  systems  were  deployed  near  the  seafloor 
by  the  University  of  Delaware  [8].  These  systems  were  ~1  km  apart  in  -100-108  m  deep  water 
(denoted  Sta05  and  Sta07  in  Fig.  1).  A  hydrophone  receiver  monitoring  the  source  transmissions  was 
suspended  from  the  R/V  Kilo  Moana  approximately  mid-way  between  the  seafloor  systems.  The 
transmissions  were  48  ms  duration,  7-13  kHz,  linear  frequency  modulated  (LFM)  chirps  that  were 
transmitted  repeatedly  every  144  ms  for  40  s.  The  two  sources  alternated  transmitting  every  2  min  for 
an  extended  period  of  time.  Examples  of  the  time-varying  channel  impulse  response  (CIR)  from  each 
source  to  the  monitor  hydrophone  are  shown  in  Figs.  1(b)  (July  10,  2200  UTC)  and  1(c)  (July  10,  2206 
UTC)  [13].  The  two  vertical  lines  represent  the  direct  and  bottom  bounce  arrivals.  The  time-varying 
arrivals  that  follow  in  pairs  correspond  to  the  surface  bounce  and  the  bottom-surface  bounce.  Note  the 
emergence  of  scattered,  slanted  later  arrivals  that  collectively  generate  a  striation  pattern  in  the  tail. 
These  patterns  run  in  opposite  directions  for  the  two  sources.  In  this  case,  the  dominant  surface  wave 
direction  was  from  the  NE  and  roughly  aligned  with  the  two  sources  and  monitor  hydrophone. 


2 


A  ray  tracing  algorithm  capable  of  treating  a  traveling  surface  wave  field  can  be  employed  to  simulate 
the  time -varying  CIR  caused  by  a  dynamic  sea  surface  [13].  Both  ray  travel  time  and  amplitude  are 
computed  as  the  point  of  intersection  with  the  sea  surface  follows  continuously  along  the  surface 
resulting  in  the  observed  striation  pattern.  A  simplified  example  consisting  of  a  sinusoidal  sea  surface 
with  ideal  half-space  below  (1500  m/s  sound  speed)  will  illustrate  the  concept  [13].  Both  the  source 
and  receiver  are  positioned  at  10  m  below  the  mean  surface  as  shown  in  Fig.  2(a).  The  horizontal 
separation  of  200  m  contains  5  full  wave  cycles.  The  simulated  time-varying  CIR  using  the  ray  tracing 
model  is  shown  in  Figs.  2(b)  and  2(c)  when  the  surface  wave  is  traveling  from  left-to-right  (forward) 
and  from  right-to-left  (backward),  respectively.  The  CIR  is  computed  every  50  ms  for  a  20  s  period. 
The  first  vertical  line  at  133.3  ms  is  the  direct  horizontal  arrival.  The  second  vertical  line  near  134  ms 
is  the  surface  bounce  from  the  vicinity  of  the  specular  reflection  point  in  the  middle  (i.e.  100  m). 
Specifically,  the  shortest  surface  bounce  path  occurs  during  the  passage  of  a  trough  near  midrange. 
Note  the  emergence  of  a  curved  striation  pattern  in  both  directions  immediately  following  the  surface 
reflected  arrival,  extended  over  1  ms  upward  (advancing  arrivals)  and  downward  (retreating  arrivals). 
The  symmetric  striation  pattern  is  due  to  the  simplified  geometry  simulated.  An  asymmetric  source- 
receiver  geometry  results  in  a  striation  pattern  that  is  dominated  by  upward  or  downward  striations 
depending  on  direction  of  propagation  of  the  surface  waves  [13]. 

The  KAMI  1  environment  now  is  simulated  [13].  A  realization  of  a  random  rough  surface  consistent 
with  a  concurrently  measured  surface  wave  spectrum  is  displayed  in  Fig.  3(a).  The  simulated  time- 
varying  CIR  at  the  receiver  is  shown  in  Figs.  3(b)  for  Sta07  source  (forward  traveling  surface  wave) 
and  Fig.  3(c)  Sta05  source  (backward  traveling  surface  wave).  The  paths  are  numbered  on  top  of  the 
figure  and  correspond  go  (1)  direct,  (2)  bottom  bounce,  (3)  surface  bounce,  and  (4)  bottom-surface 
bounce.  The  modeling  results  in  Fig.  3  can  be  compared  to  the  experimental  data  in  Fig.  1.  While  there 
are  differences  in  detail,  the  simulation  results  appear  to  capture  the  essential  features  of  the  striation 
pattern  seen  in  the  experimental  data.  Specifically,  the  striation  highlighted  in  the  box  emerges  upward 
for  the  forward  propagating  surface  wave  (moving  in  the  same  direction  as  the  source-receiver  path) 
and  downward  for  the  backward  propagating  surface  wave  (moving  in  the  opposite  direction  as  the 
source-receiver  path). 

IMPACT  /  APPLICATIONS 

Acoustic  data  communications  is  of  broad  interest  for  the  retrieval  of  environmental  data  from  in  situ 
sensors,  the  exchange  of  data  and  control  infonnation  between  AUVs  (autonomous  undersea  vehicles) 
and  other  off-board/distributed  sensing  systems  and  relay  nodes  (e.g.  surface  buoys),  and  submarine 
communications. 

RELATED  PROJECTS 

In  addition  to  other  ONR  Code  3220A  and  Code  321US  projects  investigating  various  aspects  of 
acoustic  data  communications  from  both  an  ocean  acoustics  and  signal  processing  perspective,  two 
recently  completed  MURIs  also  focused  on  acoustic  communications  (W.  Hodgkiss,  “Impact  of 
Oceanographic  Variability  on  Acoustic  Communications”  and  J.  Preisig,  “Underwater  Acoustic 
Propagation  and  Communications:  A  Coupled  Research  Program”). 
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Figure  1.  (a)  Schematic  of  the  KAMI  experiment  along  with  a  top  view  of  the  source-receiver 
geometry  and  dominant  surface  wave  direction  (NE  to  SW).  Two  sources  located  at  Sta05  and  Sta07 
were  moored  2  m  above  the  seafloor  and  a  monitoring  hydrophone  was  suspended  from  the  R/V 
Kilo  Moana  at  25  m  depth,  (b)  Channel  impulse  response  from  the  Sta07  source  (July  10,  2200 
UTC).  (c)  Channel  impulse  response  from  the  St aO  5  source  (July  10,  2206  UTC).  The  horizontal 
axis  is  relative  delay  and  the  vertical  axis  is  ping  time.  The  two  vertical  lines  represent  the  direct  and 
bottom  bounce  arrivals.  The  time-varying  arrivals  followed  in  pairs  and  correspond  to  the  surface 
bounce  and  the  bottom-surface  bounce.  Note  that  a  striation  pattern  in  (b)  and  (c)  emerges  in 

opposite  directions  in  the  tail. 
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Figure  2.  Simulation  results  for  a  sinusoidal  surface  and  a  symmetric  source-receiver  geometry  in 
an  ideal  half-space,  (a)  The  surface  wave  frequency  is  0.2  Hz  with  a  wavelength  of  40  m  and  the 
wave  height  of  2.8  m  peak-to-trough.  The  sound  speed  is  1500  m/s.  The  source  and  receiver  both  are 
at  10  m  below  the  mean  sea  surface  and  separated  by  200  m.  (b)  Modeled  channel  impulse  response 
when  the  surface  wave  is  moving  from  left-to-right  (forward),  in  the  same  direction  as  acoustic 
propagation,  (c)  Modeled  channel  impulse  response  when  the  surface  wave  is  moving  from  right-to- 
left  (backward),  in  the  opposite  direction  of  acoustic  propagation.  The  striation  pattern  emerges  in 
both  directions  (upward  and  downward)  and  is  symmetric  regardless  of  the  surface  wave  direction. 
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Figure  3.  Broadband  simulation  corresponding  to  the  KAM11  experiment  displayed  in  Fig.  1.  (a) 
Realization  of  a  smooth  surface  wave  profde  based  on  the  measured  directional  wave  spectrum.  The 
surface  wave  is  assumed  to  travel  from  right-to-left,  (b)  Modeled  channel  impulse  response  from  the 
Sta07  source,  (c)  Modeled  channel  impulse  response  from  the  Sta05  source.  Panels  (d)  and  (e)  both 
show  the  region  of  illumination  and  close-up  view  of  rays  interacting  with  the  irregular  surface, 

corresponding  to  panels  (b)  and  (c),  respectively. 
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